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ELECTROMAGNETIC WAVE PROPAGATION AND ABSORPTION IN 
MAGNETISED PLASMAS: VARIATIONAL FORMULATIONS AND DOMAIN 

DECOMPOSITION * * 

Aurore Back\ Takashi Hattori\ Simon Labrunie^, Jean-Rodolphe Roche^ and 

Pierre Bertrand^ 

Abstract. We consider a model for the propagation and absorption of electromagnetic waves (in 
the time-harmonic regime) in a magnetised plasma. We present a rigorous derivation of the model 
and several boundary conditions modelling wave injection into the plasma. Then we propose several 
variational formulations, mixed and non-mixed, and prove their well-posedness thanks to a theorem 
by Sebelin et al. Finally, we propose a non-overlapping domain decomposition framework, show its 
well-posedness and equivalence with the one-domain formulation. These results appear strongly linked 
to the spectral properties of the plasma dielectric tensor. 

Resume. Nous considerons un modele de propagation et d’absorption d’ondes electromagnetiques 
(en regime harmonique) dans un plasma magnetique. Nous presentons une justification rigoureuse 
du modele et diverses conditions aux limites modelisant I’injection de I’onde dans le plasma. Puis 
nous proposons plusieurs formulations variationnelles, mixtes ou non, et montrons qu’elles sont bien 
posees grace a un theoreme de Sebelin et al. Enfin, nous decrivons le principe d’une decomposition 
de domaine sans recouvrement, et etablissons le caractere bien pose de la formulation decomposee 
et I’equivalence avec la formulation a un seul domaine. Ces resultats paraissent intimement lies aux 
proprietes spectrales du tenseur dielectrique du plasma. 
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1. Introduction 

Electromagnetic wave propagation in plasmas, especially magnetised ones, is an enormous subject [28) . Even 
in a linear framework, the equations that describe it are generally highly anisotropic and, in many practical 
settings, highly inhomogeneous as well. The bewildering array of phenomena and parameters involved in this 
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modelling necessitates the derivation of simplified models tailored to the phenomenon under study, and to the 
theoretical or computational purpose of this study. 

Our interest lies in the numerical simulation of the propagation of electromagnetic waves near the so-called 
lower hybrid frequency in a strongly magnetised plasma. Such waves are used in tokamak technology in order 
to generate currents which stabilise or heat the plasma, thus bringing it closer to the conditions needed for 
nuclear fusion. The waves accelerate the charged particles that make up the plasma and transfer some of 
their energy to them through two main mechanisms: collisions between particles, which act as friction, and 
collisionless Landau damping. This phenomenon, caused by a resonance between electromagnetic waves and 
particles, is an efficient means of generating current in a magnetised plasma. Both mechanisms will be referred 
to as absorption. The basic physics of propagation and absorption is well understood [55]. Nevertheless, efficient 
and robust mathematical models have to be derived in order to do reliable numerical simulations in realistic 
settings. 

To perform these simulations, we have chosen to develop a finite element code which solves a suitable version 
of the time-harmonic Maxwell equations in a strongly magnetised plasma. It is thus a full-wave code in the 
plasma community parlance, as opposed to ray-tracing codes which solve the equations of geometrical optics. 
Because of their simplicity, the latter have been more popular for many years; however, it turns out that for 
most parameter regimes of practical interest, geometrical optics fails to hold |7j. This has renewed the interest 
in full-wave simulations. A full-wave code based on spectral methods in cylindrical geometry has been developed 
by Peysson et al. [53]. Nevertheless, generalising to real tokamak configurations, with an arbitrary cross section, 
requires to use the more versatile finite element method. 

On the other hand, full-wave computations in realistic settings are challenging because the lower hybrid 
wavelength is very small compared to the machine size m- This led us to incorporate domain decomposition 
capabilities in the code. This approach was preferred to, v.g., using state-of-the-art iterative methods to solve 
the huge linear system arising from discretisation, as it is known m that iterative methods perform poorly 
with strongly indefinite matrices such as those arising from time-harmonic equations. Thus, one has better split 
the computational domain into subdomains small enough to use a direct method. Another point in favour of 
domain decomposition is that the physical characteristics (such as density and temperature) typically vary over 
several orders of magnitude across a tokamak plasma. This might result in an extremely ill-conditioned linear 
system. However, this variation is normally continuous: another usual motivation of domain decomposition, 
viz., discontinuity in the equation coefficients, does not play any role here. 

In this article, we will present the theoretical and mathematical foundations of our code: the derivation of the 
physical model, a discussion of the possible variational formulations and their well-posedness, and the domain 
decomposition framework. The code itself will be presented in a future publication, with a series of numerical 
tests. A preliminary version, without domain decomposition and differing in several respects, has been reported 
in [53]. 

The outline of the article is as follows. In ([^we present a new, rigorous derivation of the model. We have felt 
this necessary, as plasma physics textbooks (such as [28] 1 generally invoke spurious assumptions, which actually 
are not satisfied in real tokamak plasmas. They start by assuming that the external magnetic field and plasma 
characteristics are (at least approximately) homogeneous, and neglecting absorption phenomena. The latter 
are only discussed as an afterthought, if at all. We shall see that inhomogeneity is a non-issue, and absorption 
can be seamlessly integrated into the model. This is fortunate, as this phenomenon plays a crucial role in the 
well-posedness of the variational formulations, with or without domain decomposition. Actually, a simplified 
model related to ours is known to be ill-posed in absence of absorption [14] , which has an extremely important 
consequence: the global heating effect does not vanish as absorption tends to zero. If the limiting model were 
well-posed, heating would be negligible when absorption is very small, as is the case in real tokamak plasmas. 
The section ends with a brief discussion of the boundary conditions that can model wave injection into the 
plasma. 
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In §33]and|3]we discuss the various possible variational formulations for the injection-propagation-absorption 
model, and prove their well-posedness. The simplest one (which we call the plain formulation) has been 
discussed in m- We recall the results of this reference, which has not been published in a journal and thus 
has not reached a wide readership. Then, we show how they extend to mixed and augmented formulations in 
the sense of m- Mixed formulations enforce the divergence condition on the electric field, and thus control 
the so-called “space charge” phenomena. Augmented formulations allow one to use the simpler nodal finite 
elements [SM ITTlfT^ etc] instead of the edge (Nedelec) elements. Furthermore, the computed solution is 
continuous (like the physical solution), which avoids spurious difficulties when coupling with other solvers 
of computational plasma physics. We conclude this part by listing a few properties of the functional spaces that 
appear in the variational formulations. 

In (J5] we present the non-overlapping domain decomposition framework for our equations. Following the 
above discussion, we focus on the mixed augmented formulation. We prove the well-posedness of the domain- 
decomposed formulation and its equivalence with the initial, one-domain formulation. This parallels and gen¬ 
eralises the work done in [3] for the time-dependent Maxwell equations in a homogeneous, isotropic, and non¬ 
absorbing medium. Notice, however, that our problem is considerably more difficult to solve numerically than 
in the latter work: the matrix of the linear system arising from the discretisation of our variational formulation 
is neither definite nor Hermitian, unlike that of [4]. 

2. Electromagnetic waves: model problem 

The physical system we are interested in is a plasma or totally ionised gas, pervaded by a strong, external, 
static magnetic field Bq{x). (We shall always denote vector quantities by boldface letters.) Such a medium can 
be described as a collection of charged particles (electrons and various species of ions) which move in vacuum 
and create electromagnetic fields which, in turn, affect their motion. Electromagnetic fields are, thus, governed 
by the usual Maxwell’s equations in vacuum: 

^ 1 95 

^~~dt’ = + ( 1 ) 

div5 = e/£o, divB = 0. (2) 

Here 5 and B denote the electric and magnetic fields; g and J' the electric charge and current densities; Eq 
and pto the electric permittivity and magnetic permeability of vacuum, with Eq = 1. 

2.1. Wave propagation equation 

The electromagnetic field is the sum of a static part and a small perturbation caused by the penetration of 
an electromagnetic wave. The latter is assumed to be time-harmonic. To simplify the discussion, we assume 
the plasma to be in mechanical and electrostatic equilibrium in the absence of the wave. Thus, the electric and 
magnetic fields can be written as 

6{t,x) = e^[E{x)e~^‘^*] and B{t,x) = Bo{x) + e^[B{x)e~^‘^*], (3) 

where i = -v/—1, 5ft denotes the real part, e <C 1, and oj > 0 is the wave frequency. In the same way, we have 

J'{t,x) = e^[J{x)e~'‘^*] and g{t,x) = e^[p{x)e~'‘^*]. (4) 

The static parts of 5, JT, g are zero by the equilibrium assumption. Furthermore, the static magnetic field 
satisfies div Bq = 0 and curl Bq = 0, as its sources are supposed to be outside the plasma. Plugging this ansatz 
in the Maxwell equations o, ©, we find: 

curl E = \uj B, curl B + iwc”^ E = po J 
dYvE = p/so, divB = 0. 


(5) 

( 6 ) 
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Eliminating the variable B between the two equations in ([5]), one finds 

curl curl E — ^E — iujfiQ J. (7) 

We will shortly show that the medium obeys a linear, inhomogeneous and anisotropic Ohm law: 

J{x) = a{x) E{x). (8) 

The expression of the conductivity tensor will be derived in the next section. Finally, Eq. © becomes 

curl curl E — f / H —E = 0, (9) 

\ eoLo ) 


where / is the identity matrix. 

2.2. The plasma response tensor 

As in [231115], the current density J in (|S]) appears as the sum of a “classical” part, which can be explained 
by a fluid model, and a kinetic correction arising from Landau damping; both are linear in E. Let us begin with 
the classical part. The particles species are labelled with the subscript the charge and mass of one particle 
are called and m,^. In a first approach the plasma is assumed to be “cold”, be., the thermal agitation of 
particles, and thus their pressure, is negligible. Each species obeys the momentum conservation equation^ 

Q'XJL 

——A + (u^ ■ V)m^ — (f + x S) + VcU^ = Q (10) 

at 

where denotes the Eulerian fluid velocity and Uc is the ion-electron collision frequency. (Collisions between 
particles of the same species do not change their bulk velocity.) The fluid velocity and the particle density 
n^(t, a:) are linked to the electric charge and current densities by: 

g = , J = := . 

Multiplying Eq. OT by we find 

^ - {J, ■ - — {g,£ + J, X B) + iy,J, = 0. (11) 

ot Qci 

Then, we use the ansatz ©-0- More specifically, for each species we assume 

J7’^(t, ai) = e 5R[Jc^(a;)e“'“’*] and g^{t,x) = q^ n^{x) + e^[p,;{x)e~^‘^*]. 

The static part of J^ vanishes, as the plasma is at rest when e = 0. At order 1 in e, one can discard the term 
altogether; in g^ £ and J^ x B, only the terms in q^ E and x Bq survive. Furthermore, we 
introduce the plasma and cyclotron frequencies for each species 


'^P‘; ■- 


I n°q‘^ 

£o 


:=—|Bo|, 
rric 


^It can be derived by integrating in velocity the Vlasov equation, see for instance m- 


( 12 ) 
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as well as = sign(g^) and the unit vector b = Thus we obtain the relationship: 

i(a; + iPc) + £ooJp^ E + S^LUa x b = 0. (13) 

At each point x, one considers an orthonormal Stix frame [5S] {ei{x),e 2 {x),e^{x) = b(x)). For any vector 
field V, one denotes ti|| = ^363 and v± = Viei + V2e2 the components of v{x) parallel and perpendicular 
to Bq{x). Taking the cross product of m with 6 on the right, we have: 

i(u; + ivc) X b + £oUJp^ E x b — 6^u!c,i J<;_l = Oj (14) 

as there holds J± = b x {J x b). Again, we take the cross product of (IT4l) with 6 on the left: 

i(a; + iuc) J^j_ + eoUjl^E± + 5 ^uJc^ J^x b = 0, 
which allows us to eliminate x b in (jTT)) : 

((w + iuc)^ - J<;± = i(w + iuc) Ejl - Exb. (15) 

The parallel current J^|| := (J,; • b) 6 is obtained by taking the dot product of m with b: 


_ i£o^p, 

'^^11 "cc + ip. '‘■ 


(16) 


Thus, the total current density J,^ = Jc;|| + J,^j^ of the species c is given as: 


U! + iUc 


^ i(a; + iPc)eoWp^ ^ 


- 


(a; + ip,)2 - ujI (c. + iue)" - ujI 


E X b. 


(17) 


Taking all species into account and setting a{x) := w + iuc(a:), we find the expression of the “classical” 
current density: 


J,., = kocc E — ^11 + - E - E E X b. 

UJCX " hi — hT^ h! — hl^ 


a 
w 


1 ^ dr iOccUjl 


uj ^— ui, 


(18) 


= :/3 


= :T 


=:<5 


In the Stix frame, we have E — Eiei + E 2 e 2 + E^b, E^^ = E^b, Exb — ^263 — Eie 2 . This gives the classical 
part of the conductivity tensor in ([5]): 


^7 —id 0 ^ 

^cia = ieow I id 7 0 I . 

0 0 /3> 


(19) 


The classical dielectric tensor is thus: 


e :=I 


SqUJ ' 


cla 


/l - 7 

id 

0 \ 


-iD 

0 \ 

-id 

1-7 

0 

:= iZl 

S 

0 

V 0 

0 

1-/3/ 

VO 

0 

p 


( 20 ) 
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where the functions S, D and P are given by 


S{x) 

a{x) ^p<;(.x) 

^ 2^aHx)-u;lix)’ 

(21) 

D{x) 

1 St^ UJc<;ix)LOp^-{x) 

UJ ^ a‘^(x) — UJe^(x) ’ 

(22) 

P{x) 

s s 

"3" I 

1 

1—1 

II 

(23) 


We proceed with the Landau damping part. As it appears [2^, only electron Landau damping in the direction 
parallel to Bq plays a significant role. The “resonant” current generated by this effect is thus of the form: 

J res(^) — Te(^) .^|| (^)- 

The coefficient 7 e is derived from a local linearisation of the Vlasov equation in the neighbourhood of the 
point X. Following the classical treatment by Landau |21] . and assuming a Maxwellian distribution function at 
order 0 in e, one finds pS] : 


7e = eow 


2 fc. \kB Tf 


rUe 


3/2 


exp - 


me 

2 kf. fee Te 


(25) 


where Te is the electron temperature (the subscript e refers to electrons), k^ the Boltzmann constant, and fey 
is the component of the wave vector fe parallel to Bq. 


Adding the two contributions, J 
matrix appearing in ([8]): 


Jc\a + Jres-i we find the expression in the Stix frame of the conductivity 

/O 0 0 \ 

5: = 571,, + 0 0 0 . 

Vo 0 7e/ 


In other words, the equation ([9]) which governs electromagnetic wave propagation and absorption in the plasma 
can be rewritten as: 

curl curl FI — FI = 0, (26) 

with the plasma response tensor given by: 

/ S{x) -iD{x) 0 \ . /O 0 0 \ 

K{x)=\ iD{x) S{x) 0 +— 0 0 0 . (27) 

V 0 0 Pix) J \o 0 je{x)J 

"- ^^ -- 

£(£c)=classical dielectric tensor Landau term 


The tensor K.{x) is not Hermitian as soon as Uc or je > 0. 

2.3. The injection-propagation-absorption model 

Let n be a bounded open domain in which represents the plasma volume in the tokamak. From the 
previous subsection, we know the propagation-absorption equation: 


curl curl E — E = 0 in 11. 

C‘^ - 


(28) 
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The entries of the plasma response tensor K_ are given as functions oi x € fl. The divergence equation 

dW{KE)=0 infl, (29) 

is a direct consequence of the previous one and may appear redundant. Nevertheless, it will play an all-important 
role in the derivation of the mixed and augmented variational formulations. 



Figure 1. A cross-section of the domain il. 

Furthermore, various boundary conditions may be considered. Let F be the boundary of fl, and n the 
outgoing unitary normal vector. This boundary is made up of two parts (see Figure [IJ: Fyi corresponds to the 
antenna and Fc = F \ F^i is the remainder. Introducing the surface current jj^ flowing through the antenna, 
the usual jump relations between media [6ll^ give: JBt = —^ where Bt denotes the component of B 
tangent to the boundary. Using the first part of ([S]), we deduce: 

curlE X n = iojfiojji on Fa. (30) 

It appears as a Neumann (natural) condition. This modelling seems more relevant than that of [2411261127] , where 
j is treated as a ficticious volumic current in U. Alternatively, one can use an essential (Dirichlet) condition: 

E X n = Ea X n on Fa, (31) 

where Ea is the electric field excited at the antenna. On the rest of the boundary, we use a perfectly conducting 
(homogeneous Dirichlet) boundary condition 


E X n = 0 onFc, 


(32) 


for the sake of simplicity. 


3. Variational formulations 


3.1. Functional setting 

In the whole article, we suppose that the domain D is a curved polyhedron, i.e., a connected Lipschitz domain 
with piecewise smooth boundary such that, near any point of its boundary, O is locally C°°-diffeomorphic to a 
neighbourhood of a boundary point of a polyhedron. This definition includes both smooth domains and straight 
polyhedra. Furthermore, the boundaries Fa, Fc are collections of smooth faces separated by smooth edges, 
possibly meeting at vertices. 
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We shall use boldface letters for the functional spaces of vector fields, v.g., L^(n) = L^(n)^. The inner 

product in L^(n) or L^(fl) will be denoted (• | ■)n- Most unknowns and test functions are complex-valued, so 

this inner product is Hermitian. “Duality” products {ip,v)v will be linear in the first variable ip and anti-linear 
in the second v; the subscript V indicates the space to which the latter belongs. In this case, ip G V, the space 
of anti-linear forms on V, which we call its dual for short. The subscripts D, V may be dropped when the 
context is clear. 

Let H(curl; D) be the usual space of square integrable vector fields with square integrable curl in D. We 
introduce the ranges of the tangential trace mapping ■ v v x n and the tangential component mapping 
TT-f '■ V 1 -^ '■= n X {v X n) from H(curl; D): 

TT(r) := {if G H-i/2(r) : g H(curl; D), ip = v x n|r}, (33) 

TC(r) := {A G : 3v G H(curl; D), A = -urir}- (34) 

These spaces have been described in [^, where they are respectively denoted Hy ' (divr,r) = TT(r) and 

H])^^^(curlr,r) = TC(r). Furthermore [TU] . they are in duality with respect to the pivot space (T) := {w G 
L^(r) : w ■ n = 0}. This allows one to derive an integration by parts formula, valid for any u, v G H(curl; D): 


{u I curlr;)n - (curlw | v)q = {u x n, r;T)Tc(r) ■ (35) 

Traces on a part of the boundary, v.g. Tc, can be defined straightforwardly: the range spaces, calledH|| ^^^(divrc, Tc) 
and H^^pg^(curlrc,Tc) in [5], will be denoted here TT(rc), TC(rc). Introducing the space 

(curl; D) := {u G H(curl; D) : m x n^-pc — 0}) (36) 


i.e., the subspace of fields satisfying the essential condition (15^ . the range of the trace mappings on the rest of 
the boundary T^ will be denoted 


TT(r^) = {c^GH-i/2(r^):3^gH^(curl;D), <^ = -cxn|r^}, (37) 

= {p G TT(r^) : the extension of by 0 to T belongs to TT(r)}, 

TC(r^) = {A G H-i/2(r^) : g H^(curl; D), A = Uyir^} (38) 

= {A G TC(r^) : the extension of A by 0 to T belongs to TC(r)}, 


instead of the “learned” notations H|j’^^^(divp^, Ta) and H(j(^^^(curlp^, T^) of [5]. The spaces TT(r^) and TC(rA) 
are in duality with respect to the pivot space (T^), and similarly for TC(r^) and TT(r^). 

Similarly, we introduce the Hilbert space: 


H(divK;D) := {u G L^{n) : div(^M) G L^)!!)}, 


(39) 


endowed with its canonical norm. If G L°°(D; Al 3 (C)), it can be alternatively characterised as H(div^; D) := 
{u : Ku G H(div;f2)}, and the conormal trace Ku • n of a field u G H(div^;D) is defined as an element 
of H~^/^(r). Then we have another useful integration by parts formula, valid for all u G H(divX;D) and 
p G H^(D): 

(div(Kti) I p)^ + (Ku I grad p)^ = (Ku ■ n, p}Hi/ 2 (r) ■ (40) 

If p G HJ(S1), the above formula can be extended to u. G L^(D). In that case, divliTw) G H“^(D) and 


(div(KM),v5)Hi(n) + (Ku I gradp)^ = 0. 


( 41 ) 
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In both cases, the scalar product (Ku \ gradt/j)^ can be written {u \ K* grad t/jlr., with K * the conjugate 
transpose of X. 

3.2. Non-mixed formulations 

Applying the Green formula (IMI) to (l28)l with the boundary conditions (130)1 and ()32l) , the electric field appears 
as solution to the following variational formulation: Find E G H^(curl; fl) such that: 

a{E, F) = l{E), V F G Hf (curl; fl) (42) 

where the forms a and I are: 

a{u,v) := (curlit I curlt!)n - ^(^M I v)n, (43) 

l{v) ■= iw/ro(j^,r>T)TC(p^). (44) 

The formulation (|42[) will be called the plain formulation. It can be regularised by adding to both sides a term 
related to the divergence. To this end, we introduce the spaces: 

X(K;II) := H(curl;n) nH(div^;II), (45) 

X%{K;n) := H^(curl;II)nH(divK;II), (46) 

endowed with their canonical norms and inner products. Using Eq. (ED), we obtain the augmented variational 
formulation: Find E G X^(^;U) such that: 

a,{E,E) = l{F), VFgX^(K;U) (47) 

with the augmented sesquilinear form Os(-, •) (s G C) defined on X(^; U) as 

as{u, v) := a{u, u) + s (div(Ku) \ div(Kt!))n. (48) 


3.3. Mixed formulations 

Alternatively, the divergence condition (1291) can be considered as a constraint. Starting from the plain 
formulation (1^ . we introduce a Lagrangian multiplier p G IlJ(r2) to dualise this constraint {cf. (1411) 1 and we 
obtain a mixed unaugmented formulation which writes: Find {E,p) G H^(curl;r2) x Hj(U) such that 

a(F,F)+^(^ = 1{F), VFGH^(curl;U), (49) 

P{E,q) = 0, VgGHj(U), (50) 

where the form /3 is defined on L^(fl) x H^(U) as: 

P{v,q) := -{Kv I gradg)^. (51) 

If we start from the augmented formulation (1471) instead, we introduce a Lagrangian multiplier p G L^(r2) 
and arrive at the mixed augmented variational formulation: Find {E,p) G X'^rlK: il) x L^(U) such that 

a.(F,F) + 6(^ = 1(F), VFgX^(K;U), (52) 

b{E,q) = 0, VqGL2(U), (53) 

with &(•,•) defined on H(div,^;U) x L^(r2) as: 


b{v,q) := (div(Ka) | q)n. 


(54) 




10 


TITLE WILL BE SET BY THE PUBLISHER 


3.4. Essential boundary conditions 

The variational formulations for the essential conditions m and (ESI) are obtained in a similar way. Using test 
fields satisfying F x n = 0 on dfl in the Green formula (1551) . one derives the plain and augmented formulations: 
Find E € H^(curl;n), satisfying (I5T]) . and such that 


a(F,F)=0, VF G Ho(curl;17); 

(55) 

Find E G X^(^;17), satisfying (|3ip. and such that 


a«(F,F) = 0, VFgX^(K;17), 

(56) 

where we have set: 


Xn{K; 17) := Ho(curl; 17) D H(divK; 17). 

(57) 


To analyse these formulations, one splits the electric field as E = Ea + E° , where Ea is a lifting of the boundary 
data and E° satishes the perfectly conducting condition on the whole boundary. Thus, one has to assume at 
least that Ea G TC(rA). This is obviously sufficient for the unaugmented formulations. For the augmented 
formulations, it is necessary to have a lifting in X(X; fl). As we shall see in Remark 14.121 the existence of such 
a lifting does not entail any supplementary condition on Ea- In this case, E° belongs to the space Xjv(^; 17). 
The plain and augmented formulations satisfied by E° are respectively: 

Find E° G Ho(curl; 17) such that 

a{E°, E) = (/, F)H„(euri;n) := VF G Ho(curl; 17); (58) 

Find E° G 'K.n(K\ 17) such that 


as{E\F)=Ls{E), VFgX^(K;17). (59) 

The mixed augmented formulations satisfied by E and E° write: 

Find {E,p) G X^(^;17) x L^(17), with E satisfying (1511) . and such that 

a«(F,F) + 6(^ = 0, VFgX^(K;17), (60) 

b{E,q) = 0, V9GL2(17). (61) 

Find {E°,p) G XvliT: 171 x L^(17) such that 

as{E°,F)+^F;^ = Ls{F), VFgX^(K;17), (62) 

b{E\q) = £(g), VgGL2(17). (63) 


The right-hand sides are given by: 

Ls{v) := -a{EA, v) - s (dWiKEA) \ div(Ku))n := (/, ti) -|- s (5 | div{JKv))n ; (64) 

£{q) := -{div{KEA) \ g)n := {g \ q)n ■ (65) 

The reader may write the mixed unaugmented formulations as an exercise. 

As usual when dealing with non-homogeneous essential conditions, we shall use the formulations in E°, 
Eqs. dSll), dSH), and (I55|) - (I551) to prove well-posedness. However, in practice, we discretise the formulations 
in E, Eqs. dnn), dini), and (EOll-dnil): both conditions m and (1551) are handled by a pseudo-elimination 
procedure following a local change of basis [18] . 
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4. Well-posedness of the problem 

In this section, we summarise the results of m — which deals with the plain formulation when 7 e = 0 — 
and show how they extend to our various formulations. We shall make the following assumption throughout 
the article. 

Hypothesis 1. The real functions Vc, ajc<; and ojp^, for each species (ions and electrons) are bounded above 
and below by strictly positive constants on O. The function 7 e is non-negative and bounded above. 

Remark 4.1. The collision frequency Uc is given by the following expression m, where Z is the ion charge 
number (i.e., their charge is equal to Z \qe\): 


Vr=\ — 


'^pe 


In A 


A 


with: A = 


1271 


Eq ^B Te 
Ueql 


3/2 


( 66 ) 


A plasma is characterised by A 1. In this framework, and recalling the expressions (ED of UJcc;^ ojpc; and (ESD 
of 7 e, one checks that Hypothesis [T] is satisfied provided the densities and the electron temperature Tg are 
bounded above and below by strictly positive constants on H. This is the case in all practical settings. 


4.1. Spectral properties of the plasma response tensor 

It is not difficult to check that the eigenvalues of the matrix R are 


Ai — S D, X 2 — S — D, A 3 — P -\ - 7 e’ 

SqUJ 

Furthermore, K_is a. normal matrix {Kf K_ — K K *). and its singular values are the moduli of its eigenvalues. 
Then, one deduces from dUD-dMD the expression of the imaginary parts (SAi)i=i... 3 : 


SAi(a:) 

5A2(ai) 

3A3(a;) 


J/c(£c) 


UJ 


E 




«(ai) - + u 2 ( a .))2 + 4 ^ 2 ^ 2 (^) 


Uc(®) 

^ - a ;2 + u 2 ( a ;))2 + 


uj^,{x) 


UJ 


<; 

Vc{x) 


^cc2 + u2(ai)) 




[(w - 5^u:c<i{x)f 
[(w + 5<;UJc<;{x)Y 




From the above calculations, one easily infers a fundamental bound. 

Lemma 4.2. Under Hypothesis\^ there exist two constants q > C > 0, dependent on oj, such that 

q{z*z) > \z*K_{x)z\ > 3 [( 2 :*R(a:) 2 :)] > ((z*z), Mz £ C^, Mx G H. (67) 


Remark 4.3. When the ujcq and LOp^, Vc and 7 e have typical values for tokamak plasmas, and w is of the 
order of the lower hybrid frequency, one has 5RAi > 0, while KA 2 < 0 and KA 3 < 0. No lower bound holds for 
mz*K{x)z)]\. 

4.2. Coercivity and inf-sup condition 

We recall the fundamental result of Sebelin et al. Ell- 

Theorem 4.4. Let V and H be Hilbert spaces such that the embedding V ^ H is continuous. Let a(-,-) be a 
sesquilinear form on V x V. If there exists three strictly positive constants a. A, 7 such that: 
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(1) the real part o/a(-,-) is Carding-elliptic on V, i.e.; 

|K[a(z;,p)]| > a||v|l^ - AIIpIII^, Vw e F ( 68 ) 

( 2 ) the imaginary part o/a(-,-) is H-coercive, i.e.; 

|3[a(p,'y)]| > 7 ||z;||^, Vv € V, (69) 

then the sesquilinear form a is V-elliptic. 

Combined with Lemma SSI this theorem shows the well-posedness of the non-mixed formulations, by Lax- 
Milgram’s lemma. 

Theorem 4.5. There exists a unique solution to the plain formulations (1421) and dSll), hence to (i55|) . The 
same holds for the augmented formulations dUl) and dSH) — and thus for (EH) — provided > 0 and 3s < 0. 

Proof. As in one uses Eq. (l67l) to check that the form a given by (l43|) is continuous onV = H(curl; fl) and 
satisfies the assumptions of Theorem 14.41 with H = L^(n). Thus, it is coercive (and continuous) on H(curl; 17), 
a fortiori on the closed subspaces Ho(curl;17) and (curl; 17). When Ks > 0 and 3s < 0, the same applies 
to the form Og given by (H51) on V = 'K(K-. 17). This form is coercive and continuous on the closed subspaces 
of X{K; 17), v.g., '^n{K; and X.%{K; 17). □ 

To prove the well-posedness of mixed formulations, we have to check an inf-sup condition. This can be done 
by following the lines of Enna. 

Proposition 4.6. The sesquilinear form /3 defined by (ED) satisfies an inf-sup condition on Ho(curl; 17) x Hi(17), 
i.e.there exists Cg > 0 such that 

VqGHj(17), sup > C/j Ikllni • (70) 

vEHo(curl;r2) ||”^|| H(curl) 

Proof. Fix q G Hj(17) and set v — gradq G Hq (curl; 17). Lemma 14.21 shows that 

\0{v,q)\ = \{Kv I gradg)| = \{Kv \ t!)| > CII'^IIl^ = CII'^IIlHI gi’adg||L 2 . 

1/2 

On the other hand, ||v||H(curi) = (II^’IIl^ + II curluH^^) = ||t;||L 2 , and || grad( 7 ||L 2 > CHgllHi by Poincare’s 
inequality. Hence the conclusion. □ 

To proceed to the mixed augmented case, we state and prove a useful lemma. 

Lemma 4.7. For any f G H~^(17), the elliptic problem: Find (p G Hj(17) such that 

- Ax()) :=-div(Kgrad())) = / (71) 

admits a unique solution, which satisfies I'i/'Ihi ^ C'II/IIh-i /a*' some constant C. 

Proof. Using (|411) with u = grad p, the variational formulation of (l711) writes: 

a(((), ■!/;):= (Kgrad 0 I grad■!/;)n = (/,'0)Hi(n), Vi/; G HJ(17). (72) 

By Eq. (1571) . the form a satisfies 

’llV’Inqo) > > CIV'lHi(n)> e Hj(17), 

i.e., it is continuous and coercive on H()(17), and the formulation is well-posed by Lax-Milgram’s lemma. □ 



TITLE WILL BE SET BY THE PUBLISHER 


13 


Remark 4.8. Elliptic problems with the operator A^* are also well-posed. 

Proposition 4.9. The sesquilinear form b defined by (1541) satisfies an inf-sup condition on 'X.n(K', fl) x L^(n), 
i.e.there exists Cb > 0 such that 


Vg G L2(11), 


sup 

vGXiv(i£;n) 


ll^llx 


>aikiiL- 


(73) 


Proof. Fix q G L^(0). According to Lemma IdTZl there exists (j) G Hj(r2) such that A^()) = q. Setting v = grad (f, 
we have v G Ho(curl; fl) and divlJCul = q, hence v G H(div^; ft) and finally v G Xjv(^; H). It is bounded 
as: 


ll^’llx = II^’IIl^ + II curlu||^2-h II divKu||^2 

= ||grad(/)||L2-tO-h||divKgrad((i||L2 

= klH^ + lkllL^ < (i + c^^)lkllL^, 


On the other hand 
Finally 


which we had to prove. 


IK^^,9)I = \{divKv I 9)1 = 1(9 | 9 )] = || 9 || 2 , 


l^(^^g)l > 

ll^llx - 


VTT^ 


Il 2 


□ 


Theorem 4.10. There exists a unique solution to the mixed unaugmented formulation (|4^-(|5nil. and to its 
counterpart for the Dirichlet boundary condition. The same holds for the mixed augmented formulations (j52|l - 
dsai) and (|62l)-(p|l, provided 5Rs > 0 and 3s < 0. Thus, the problem (|60|-(|6T1l is well-posed in this case. 


Proof. The forms a and Og are coercive, in particular, on the kernels of the forms /3 and b. The form b is 
obviously continuous, and so is /?(•, •) thanks to the boundedness of the entries of X. The inf-sup conditions 
dlol) and (l7^ are exactly those needed for the Dirichlet formulations. In the Neumann case, they remain valid 
when replacing Ho(curl;n) with H^(curl;n) or XAr(^;r2) with X.'trlK: Dl. as the supremum is greater on 
the bigger space. We conclude by the Babuska-Brezzi theorem. □ 


To conclude this subsection, we observe that all formulations are equivalent to one another. For instance, 
the unique solution to the plain formulation satisfies (j28|) in hence div KE = 0 and E G 'SSfiiK: Dl is 

solution to the augmented formulation. Similarly, {E, 0) is solution to both mixed formulations, thus it coincides 
with their respective unique solutions. 


4.3. Miscellaneous properties 

Here we collect and discuss some useful properties of our functional spaces. First, one has a Helmholtz 
decomposition of vector fields into gradient and “j^-solenoidal” parts. 

Lemma 4.11. For any u G L^(r2) there exists a unique pair {(j),UT) G H(j(r2) x L^(f2) satisfying the conditions 


u = grad (j) -\-ut, div(KuT ) = 0 : (74) 

II grad((i||L2 < C||w||l2 , < C'||m||l2 . (75) 

Proof. If a solution exists, then A^(f) = divlJCn): the latter function belongs to H“^(r2) under Hypothesis [U 
with II divfiTwlllTT-i < C ||m||l 2 . Lemma 14771 shows the existence and uniqueness of (p and ut = u — gradi)), as 
well as the bounds (1751) . □ 

Remark 4.12. Obviously, gradi^ G Ho(curl;H) and ut G H(div^;r2). As particular cases: 
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(1) If It G H(curl;r2), then ut G H(curl;r2) and thus ut G Furthermore, ut x n = u x n: 

the ranges of the mappings 7 t and Try from X(^; fl) are identical to the ranges from H(curl; 17), i.e., 
TT(r) and TC(r). 

(2) If M G X(X; 17), then gradt); G Xjv(:^; 17), and the decomposition (1741) is continuous in X norm. 

(3) As a consequence of the two previous points, both Ut and gradt/) belong to Xjv(^; 17) if u does. 

Remark 4.13. Thanks to Remark 14.81 one also has a decomposition into R*-gradient and solenoidal parts: 
for any u G L^(17), there is a unique pair (()), tty) G Hj(17) x L^(17) such that 

u = K * grad (j) + ut and div ut = 0. (76) 

The above results allow one to prove two powerful theorems on the space X7v(R;17). They parallel the 
well-known results valid for scalar or Hermitian positive definite dielectric tensors. The proofs are similar to 
these classical cases and can be found in [TB], so we will not detail them here. 

Theorem 4.14. 7/17 is Lipschitz, the space Xjv(R;17) is compactly embedded into L^(17). 

Proof. Follow the lines of Weber [SD] , using Lemmas 14.2114.71 and 14.111 □ 

Remark 4.15. If Uc = 7e = 0, the proof breaks down: without absorption, one cannot establish a Fredholm 
alternative for the model of 1)2.31 with Dirichlet boundary conditions, see also Remark 14.31 With Neumann 
boundary conditions, the embedding X^(R;I7) L^(f7) is not compact when F^ ^ 0 , whatever the matrix 

field R. Thus, all usual strategies for proving well-posedness fail in the absence of absorption. Actually, there 
is every reason to believe that the model is ill-posed in this case (see Ell- 

Theorem 4.16. Assume that 17 has a boundary, and that the functions Vc, le, and (for each 
species <;) belong to C^(I7). The space Xjv(R;I7) is algebraically and topologically included in H^(f7). 

Proof Following the lines of Birman-Solomyak [^, one shows that any u G 'Kn {!£.', f7) admits a decomposition 

u — ubs + with: mbs G Hj!,(I7), G Hj(f7), G L^(I7). 

The usual elliptic theory [12], valid for the operator —A_i^ thanks to Lemma l4^ shows that (p G H^(f7), given 
the smoothness of ft and Hence u G H^(f7). 

In other words, there holds XAr(R;I7) C H)Y(f7) := {to G H^(f7) : w x n|r = 0}; the converse inclusion 
is obvious as R G C^(I7; Al 3 (C)). Furthermore the embedding H)y(I 7) ^ Xjv(R;I7) is continuous; thus the 
converse embedding Xjv(R; f7) ^ H)y(I 7) is continuous by the open mapping theorem. □ 

Remark 4.17. Under the hypotheses of the above theorem, it is thus possible to discretise straightforwardly 
the augmented and mixed augmented variational formulations of El with nodal (Lagrange or Taylor-Hood) 
finite elements. Note that this does not apply when the boundary is not smooth and has re-entrant corners, 
due to the singularity of the solution m- 

5. Non-overlapping domain decomposition framework 

For the sake of simplicity, we assume from now essential boundary conditions, and we consider (c/. 1)3.41) the 
following model problem: 


curl curl FI — ^KE 

= f 

in 17, 

(77) 

div(RF;) 

= g 

in 17, 

(78) 

E X n 

= 0 

on F, 

(79) 




TITLE WILL BE SET BY THE PUBLISHER 


15 


where the data {f,g) satisfy the compatibility condition div/ = —^ g. The mixed augmented variational 
formulation reads: 

Find (E^p) G XAr(^;f2) x L^(n) such that 


as{E,F)+b{F,p) = L,(F) := (/ I F)n + s (5 I div(KF))o, VFGX^(K;fl), (80) 

b{E,q) = ^{q)■.= {g\q)n, Vg G (81) 

As shown by the above notations, we have assumed f G L^(r2) and g G L^(n). According to dH this problem 
admits a unique solution {E,p) G X.n{K_; fl) x L^(fl), with p = 0. 

5.1. Strong formulation 

We introduce a non-overlapping domain decomposition 


_ Na _ 

n = U n, ; fliCn, i = l,...,Nd; fl* n = 0 (82) 


The exterior boundaries of subdomains are denoted Ti = T fl dfli, i = 1,..., Nd, and the interfaces between 
them Yiij = dfli fl dflj. We shall write i A j whenever is a non-empty topological surface, i.e., it has a 
non-zero area. To keep things simple, we assume that the T^ and Yij are smooth when they are topological 
surfaces, and then that the T^ n Yij and n Yi^k are smooth curves. This is generally achieved in practice. 
However, the skeleton of interfaces S = IJ^ Yij is not smooth, as there generally are (curved) dihedral angles 
between interfaces. 

The principle of domain decomposition for Maxwell’s equations has been known for some time, both in the 
time-harmonic [nimiiiiiaiis] and time-dependent 011 ] versions. Consider the solution E to (I77|)-(I7n|), and 
set Ei := E^q.. Clearly, each Ei satisfies: 


curl curl 


. - ^ke. 

= f^ in Hi, 

(83) 

d\Y{K^Ei) 

= g^ in rti, 

(84) 

Ei X n 

= 0 on Tj, 

(85) 


where fi,gi) are the restrictions of ( K . f, g) to 12^. In addition, we have the following interface conditions. 
As E G H(curl; H) satisfies (l77l) in the sense of T>'{Q), there holds: 

Ei X rii = —Ej X Tij on Yij , 

curl Ei X rii = — curl Ej x rij on Yij , 

where rii is the outgoing unit normal vector to dfli. Similarly, the condition E G H(div^, H) or Eq. (l7^ imply 

K.jEi ■ rii = —K_jEj ■ rij on Yij . 

Denoting as usual j = Vi — Vj (where i is the larger index) the jump of v across Yij^ the above interface 
conditions can be rewritten in the following way: 

[E X njs. ^ = 0, [KE ■ njs,,, = 0, (86) 

[curlEl X = 0. (87) 

Conversely, if the vector fields (£'i)i=i Nd defined on satisfy Eqs. (|55)) " (I571) in the suitable sense, the field 
E dehned on D by glueing them is solution to (IZZI)-®. 
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Proposition 5.1. Assume that the entries of K_ are continuous on Q,. As in the interface conditions (1861) 
are equivalent to = 0. 

Proof. The first condition implies = Xij rii for some scalar field defined on Denoting the 

value of ^ on this interface, the second part of (1861) then gives Xi^j [n^ • (' iiT ,- = 0. As rii is a real vector, 

Eq. (1571) then implies Xij = 0, i.e., [E'lsij = 0. The converse implication is obvious. □ 

5.2. Variational formulation 

Let us now introduce a variational formulation for the multi-domain equations (|MD“(|H71)- The mathemat¬ 
ical framework of domain decomposition for unaugmented Maxwell formulations is classical mill]. Roughly 
speaking, a vector field Ui G H(curl;f2i) iff it admits an extension u G H(curl;f2): in this respect, H(curl) 
spaces behave like the usual Sobolev spaces. The case is less straightforward with augmented formulations, even 
when ^ = J as in [4]- A field Ui G H(curl, div; 17^) does not necessarily admit an extension in H(curl, div; D); 
if it does, it is actually of regularity, at least away from T^ when this boundary is not empty. A similar 
phenomenon occurs in our case. As said in the introduction, we shall focus on the mixed augmented formulation. 

We consider the following functional spaces associated to the domain decomposition (1571) . They are endowed 
with their canonical “broken” norms. Conditions on the exterior boundary T^ are void if Tj = 0. 

Vo = {v G L^(fl) : Vz, Vi := v\q. G H(curl; Di) and Ui x n = 0 on Ti}, ( 88 ) 

W)v = {z;iGH(curl;D,)nH(div^;D,):t;iXn = Oonrj, (89) 

Wn = {z;GL 2(D) :Vz, tiio, G WV}. (90) 

Let E G XAr(^;D) be the solution to (1771) - (l79)) . and (£'i)j^i its decomposed version. Obviously, Ei G 
H(curl;Di) (~l H(divX;r2i), and it satisfies (I551) - (l55l) as argued above. Applying the Green formula (1551) 
on each subdomain and using the first-order interface condition (I87L we obtain the following formulation of 
Problem (p|) - ()87|) : 


Y,a^,s(E^,E,)+h{F,,p,) = ^L,,,(F,), VFgX^(K;D), (91) 

i i 

Y,h,{E,,q,) = ^£,(ft), V(7GL2(D), (92) 

i i 

[E X n]E., = 0, [KE ■ njs., =0, Vz A j. (93) 

The domain-wise anti-linear and sesquilinear forms Hi are defined as: 

a^^.s{u^,Vi) := (curlwi | curluj)n, -|- s (diY(Kui) \ diy{Kvi))Q. - ^(Ku^ \ Vi)Q. , (94) 

b{vi,qi) := {div{^Vi) \ qi)n, , (95) 

Li,s{vi) := (/j I Vi)n. -f s (51 | div{Kvi))ni , (96) 

■= {9i I g*)ni ■ (97) 


In order to dualise the zeroth-order interface conditions (155)) . we introduce various spaces of traces and jumps. 
As a first step, let: 


:= G H-i/2(e) : 3z; g Vo, = b X njs}. 

The notation [tixnjs stands for the ordered collection of jumps {[d x ^ Each jump belongs to TT(Eij) 

but in addition they have to satisfy some compatibility conditions j9]. This motivates the following definition. 
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Definition 5.2. The space TT(Ei_j) is made of the fields i~p^ ^ G TT(Eij) such that their extension by 0 
to S is the trace of a field in Ho(curl; fl): (f = v x 

Lemma 5.3. There holds: 

0 TT(E,,,)CS^. 

iAj 

Proof. Choose any interface and ^ S TT(Ei_j). It can be lifted to a field Vi G H(curl;ni) such that 
Vi X Tii = ipi j on Eij and x = 0 on dTli \ Tiij. Setting v = Vi on ili and 0 elsewhere, we have 
V G Ho(curl; 17) and = [t; x njs. Repeating the process for all interfaces yields the conclusion. □ 

Then, one defines the space C x as the range of the jump mapping: 


w I—s> [M]e := {[K.W ■ n]E, [tu x njs) • 


As in the case of , the jump on E is defined by the collection of conormal and tangential jumps on the 
which have to satisfy some compatibility conditions. Under the assumptions of Proposition 15.11 the data of 
jumps in this form is equivalent to that of three-dimensional ones. 

One is led to introduce a new Lagrangian multiplier A G (S|^)', and we obtain the following variational 
formulation: Find {E,p,\) G W^r x L^(r7) x (S§^)' such that 

+ yPGWM, (98) 

i i 

J2b^iE^,q^) = Y.AFi), Vge 1^(17), (99) 

i i 

(/X, [[£;]]e)sw = 0, V/XG(S^)'. (100) 

The duality products between and its dual can be expressed as 


(A, [[F]]s)§w 


(A„,[KF.n]s) + (AT,[Fx n]s) 

E , \F X n]s„,) 


iAj 

Ni 


E! ■''^*)H-i/2(aOi) F {XT,Fi X '>T'i)TT(ani) 


( 101 ) 


On the first two lines, the dualities hold between the suitable spaces. Furthermore, on any interface the 
sum of the contributions of f7i and 17^ amounts to a jump, as the normals have opposite orientation: rii = —rij] 
hence the third line, where by convention A„ = 0 on T^. 

Remark 5.4. Under the assumptions of Proposition EH the interface condition (jlOOl) is equivalent to: 

(m, [-E]e)s-= 0 V/xG(Sr)', ( 102 ) 

where S^, the space of three-dimensional jumps of fields in W^v, is isomorphic to If the matrix ^ G 
C^(f7;Al 3 (C)) and dfl is of regularity, then XAr(^;I7) C H^(I7) by Theorem 14.161 and the three- 
dimensional jump is defined in H^/^(E). 
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Remark 5.5. The multi-domain variational formulation (ra-dTuni) leads to a non-overlapping domain de¬ 
composition method. If we use nodal Taylor-Hood finite elements to discretise (l^-(IT00l). as discussed in 
Remark 14.171 we obtain a saddle-point-like linear system where the unknowns are the nodal values of {E,p, A). 
Mimicking Gauss factorisation, we obtain a new linear system, a generalised Schur complement system, where 
the unknowns are the nodal values of the Lagrange multiplier A only. To solve this new non-Hermitian re¬ 
duced system, we use a preconditioned GMRES iterative method. This algorithm induces at each iteration the 
resolution of a linear system corresponding to the discretisation of the variational formulation in each subdo¬ 
main, as in [4]. Then we obtain a non-overlapping domain decomposition method at the discretised level. A 
preconditioned direct method is used to solve the linear system in each subdomain. 


5.3. Well-posedness 

We now prove directly the well-posedness of the decomposed variational formulation. 

Theorem 5.6. The decomposed formulation (19811 - (I100|l is well-posed in Wat x L^(r2) x thus it admits 

a unique solution. 

Proof. The equations (j98p - (ll00l) can be written in the form of a mixed problem: 

Find {E,p,X) e Wjv x L2(n) x (S|f)' such that 

As{E,E)+B{E-p, A) = Cs{F), VEgWn, 

B{E-q,pi) = £(g), V(g,p)eL2(fl)x(Sg')', 

with 


As{u,v) 

B{v-q,n) 


i 

J2biivi,qi) + {fj-, [N]E)sg'- 

i 


Let V G Wat and (L’i)i=i,...,Ard be its decomposed version. Applying Theorem 14.51 in each one finds 
ai^e{vi,Vi) > Vi llvill^i ; thus Ae{v,v) > (min^ Vi) This holds in particular for v G kerS. 

We denote S = To prove an inf-sup condition, we choose {q, n) G L^(r2) x S' and seek v G Wn such 
that 


B{v,q,fi) >Ci 3 \\v\\w^ (IkllL^ + (103) 

with Cg independent of q and fi. To begin with, S is a space of traces, so its canonical norm is: 

||(p||s = inf {||i(i||w„ : [Kw ■ n]^ = pn and x njs = . 

The definition of the dual norm writes: 


||/x||§/ = sup 


KM,y)sl 

llv^lls 


(104) 


We introduce the decomposition S = ker fi © Ct^o, where verifies (/r, <Pq) = 1 and T kerp. Using (11041) . 
we deduce ||<Polls = 

Consider the continuous anti-linear form 1^ on Wat defined as 




(m> [M]s)s ; 


(105) 
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obviously, it satisfies ||^^||w^ < IImIIs'- On the other hand, a standard argument shows the existence of a 
continuous lifting from S to Wjv : 

[Rip X n] = (p-f-, [^Rip-n]=pn and ||i?¥ 3 ||w„ < O/? HtpHs. 

Then, we introduce the decomposition Wjv = ker(^^) © Cwq, with Wq = a^Rip^ and ag G C. The element 
Wq G W^v is normalised by the condition = ||/x||g,. This gives (/x, agi^o) = IImIIso hence ag = ||/x||g/, 

and finally ||i(;g||wN < Cr ||/x||s'. 

Next, consider ^ G Hj(n) solution to 

^k4> = / S L^(n), with fi = qi — divlJCwg,! in fli. 

This function is bounded as: 

l</'lHi(n) < Cl (lkllL2(n) + ||mo||w„) < C[ (||q||L2(n) + IImIIs')- 

The vector field v := Wq + grade/) G Wat satifies div Kvj = qi in Hi, and is bounded as: 

ll^llwjv = Ho+ grad</)||^2(n) + + lkllL2(n,)] 

i 

< 2[||’u;g||^^ + |(/)|Hi(n)] + lkllL2(n) ^ 02(11911^2(0) + ||/x|||/). 

But grad </> G Xjv(^; fl), which implies: 

\K v ■ n]s = \Kwn ■ n]^. and [d x n]^ = [lOg x n]^. 


We conclude that 


Biv, q, /x) = '^{divKv^ \ qi) + (/x, [MJs/s, 

i 

= ^ lki|lL 2 (Oi) + ' (/X, [[•U;o]]s)s , 

<©,'i«o>w„ = IImIIs/ 

= lkllL 2 (n) + Hlli' 

1 „ „ 


> 






N iiiyiiL2(o) 


1/2 


The well-posedness of the formulation (IMll - (llOO|) follows from the Babuska-Brezzi theorem. □ 

In order to interpret the decomposed formulation (155)) - (11001) . we shall need the following lemma. 

Lemma 5.7. A continuous anti-linear functional Lyy on Wat vanishes on 'K.n{K\ fl) if, and only if, it is of 
the form (11051) for some /x G (S§^)'. More specifically, there exists a unique pair (/x„,/x-|-) G x TC(S) 


such that: 


Lw{F) = X! / + Mt ■ [-F X rxjs.^} dcr. 


(106) 


iAj ' 


Proof. Let L\y G then Hahn-Banach and Riesz’s theorems {cf [8], Thm VIII.13) show that there exist 

Qq G L^(r2i), g\ G L^(r2i) and such that: 


VF, gW)v, L\^{F)= {gi,-F +g\ div KF,+gi,-curl F,) dH. 
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Since Wat = 0j W]y, we have (Wjv)' = It follows that any anti-linear form on Wjv can be written 


Lw{F)=^ / [qq-F + gid\YK_F + g^- cuyXF) dil, 

i=l 


with Qq G L^(n), gi G L^(fl) and ^2 ^ L^(n). We perform a Helmholtz decomposition of gg in ■K'*-gradient 
and solenoidal parts (Remark l4.13l) : 


gg = ^* gradi/t + gj^ := gj^ + gj^, with i/) S Hq(H) and divgy=0. 


Assume that Lw vanishes on Xjv(^;H). Let F G T>{Vl), and consider its Helmholtz decomposition into 
gradient and X-solenoidal parts: 

F = gvadcj)F t ■= FL + Ft, with (/> S Hq(H) and div-RTF^ = 0. 

Using OT, one immediately checks that {g^ \ Ft)q — 0 and [gT \ Fl)q = 0. Furthermore, Ft and Ft belong 
to X/vlRT: H), by Remark 14. 121 Since Lw vanishes on H), we deduce: 


0 = Lw{Fl) = 


g^ ■ Fl + gT -Fl + gidiM K_Fl 4- g 2 ■ curl F t ) dH. 

=0 =0 


By adding • Ft + gi div KFt) dH = 0, we have 



■ F + gi div KF) 


dn = 0 , 


VF e D(H). 


(107) 


This yields ,^*gradgi = in 'F>'{n). As g^ G L^(f2) and (K*) ^ G L°°(H; Al 3 (C)) by Lemma 14.21 we 
infer gi G Furthermore, gradgi = gradi/j in fl; as fl is connected, this gives gi = i/i -|- Ci, for some 

constant Ci. In particular, gi = Ci on the boundary F. 

Let w G Hj(fl) such that G L^(fl); then gradw G fl), and we have 


0 = L]v(gradw) 


Hoi 



• grad w + gi Ai£w) 


dfl 


(Sl ■ gradw 


gradgi ■ gradw) dfl 

V ^ 

=0 


+{gi,Kgradw n)H-i/ 2 (r) 


Taking w such that (l,Xgradw • n)H-i/ 2 (p) ^ 0, one deduces Ci = gi|p = 0, i.e., gi G Hj(fl). On the other 
hand, we have 


0 = Lw{Ft) 


gp • Ft -b gT ■ Ft + gi div KFt + g2 • curl Ft 

=0 =0 


dfl. 


By adding {gT ■ Fp -b g 2 ■ curl Fp) dfl = 0, we have 



• F -b g 2 • curl F) 


dfl = 0, 


VF G X>(fl). 


(108) 
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Thus, we have curlg 2 = —Qt T>'{n) and therefore then ^2 G H(curl; fl). Adding (11071) and (11081) . the form 
Lw is equal to: 


Nd . _ _ _ _ 

Lw{F) = / ((/i div KF + K* gradgi ■ F + g 2 - curli^ — curlg 2 • F) dfl. 

i=l 

with gi S HQ(n) and ^2 £ H(curl; fl). Using the Green formulas (HOI) and (1551) in each we deduce 


Nd „ 

Lw{F) = X! / (51 (SF~^ + Sf 2 T • {F X n)) dm 

i=i Jan, 

Each integral is understood as a sum of two duality products: the first between and 

the second between TC(5ni) and TT(9r2i). On an exterior boundary T^, one has = 0 and F x n — on. 
an interface the sum of the contributions of Hi and Q.j amounts to a jump, as noted in (IIOII) . Finally, we 
arrive at (I106L where is the trace of gi on E, and /x-p € TC(E) is the tangential component of 

^2 on E. These characterisations allow one to consider their restriction to each interface: and 

G TC(Eij) on each Sij. Of course, restrictions to neighbouring interfaces satisfy suitable compatibility 
conditions. 

To prove uniqueness, it is enough to show that Lw{F) = 0, VF S Wat, implies = 0 and /x-p = 0. First, 
take g € H“^/^(E), and introduce 0 € 11^(17) solution to the following variational formulation, with the form a 
from ((7^ : 

'ip) = {g, V'|e)hi/ 2 (e) G Hj(f7), 

which is well-posed as in Lemma |4.7I Performing an integration by parts in each and adding as before, we 
see that (p satisfies: 


—Ak^P = 0 in each , [(p]^. = 0 and [^grad (p ■ rx]^ = 5 on E. 

Setting F = grad (p, we have F G Wat, \KF ■ n]s = g and [F x rx]i; = 0. So: 

0 = Ij\\r{F^ = j liji [ F' F ■ ?xj ^ dfJ = (/Xni . 

As g is arbitrary, one deduces ^„ = 0 in H^/^(E). In particular, taking g supported on one interface E^^ , one 
finds = 0 in H^/^(Eij). 

For the tangential part, take G S^- By dehnition, there exists v G Vq such that [xi x n]x; = <p. In each 
subdomain introduce 

(pi G Hj!)(ni) solution to: — Aj^cpi = div(Kvi) G and Fi = Vi + grad (pi. 

There holds Fi G H(curl;r2i) D H(div^,Up) and Ft x rii = Vi x rii on 90^. Therefore, the global field 
F = {Fi}i^n.. satisfies F G Wat and [F x n]^ = [t; x njs = ip. As Lw{F) = 0, this implies 

(Mt><p)s^= 0 , V(pGS|j, be., /x-p =0 in (S|j)'. 

In particular, taking (p G TT(Eij), its extension by 0 to E belongs to by Lemma [5T^ and we infer /jPp' = 0 

in TT(E,jy = TC(E,,,)- □ 
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Theorem 5.8. The decomposed formulation (I98I) - (I100I) and the original mixed augmented formulation (1801) ^ 
(l8Tl) are equivalent: {E,p,\) is solution to (IM1) - (I100I) iff{E,p) is solution to (l80l) -dSl ]) . and 

Xn = 0, At = (curl£;)T|E • (109) 

Proof. Let {E,p,X) be the solution to (I^ - (I100I) . From (llOOl) . we have the jump conditions (1551) . and E € 
XAr(X;r2). Taking a test function F G XAr(X;n), the term (A, [[F]]i;)gw vanishes in (IMl) . which gives (l80l) . 
Then (|99p is identical to (I5T|) . This means that {E,p) G Xl) x L^(n) coincides with the unique solution 

to (I5U1)-(I5T]). 

Conversely, let {E,p) be the solution to (1501) - ((STl) . As F G X/yfFT: fl), we have automatically [F x n]s = 0 
and [KE • n]s = 0, which implies (IIOOII . As for (l99ll . we have 

^(divF, I Qi) = (divF \q) = {g\q) = J^igi \ qi). 
i i 

Define the continuous anti-linear form Lw on W: 

Lw '■ F I > 'y ) {—ai^s[Ei, Fi) — bi{Fi,pi) -\- Li(Fif ), 


which vanishes on XAr(^; D). By Lemma 15.71 there exists a unique A G )' such that 

Lw{F) = {A„ [KF-nJs + Ay • [FIT^Ie} da. 

So, Eq. (1051) is verified. On the other hand, we have remarked that the solution to (1501l - (15T1) satisfies divfjTFl = 
g and p = 0; thus the strong form of (1501) becomes: 

curl curl F-^^F = / in X*'(D). 

As a consequence, curlF G H(curl;r2). Starting again from (1051) . using the Green formulas (l35l) . (1401) in 
each fli, and taking the above equalities into account, one obtains: 

(- curl F, [F X n]E)sv -h (A„, [KF ■ n]E)H-i/2(E) + ('^t, [F x n]^)sv = 0. 

Thus we get the expressions of A„ and At- □ 

Acknowledgement: The authors wish to thank the anonymous referees for their useful remarks and suggestions. 
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